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Let A" be a closed convex polyhedron defined by a finite number of linear in- 
equalities. In this paper we refine the theory of abstract tubes (Naiman and 
Wynn, 1997) associated with K when K is perturbed. In particular, we fo- 
cus on the perturbation that is lexicographic and in an outer direction. An 
algorithm for constructing the abstract tube by means of linear programming 
and its implementation are discussed. Using the abstract tube for perturbed K 
combined with the recursive integration technique proposed by Miwa, Hayter 
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and the distribution functions of studentized range statistics are exhibited as 
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1. Introduction 

Let A = (fli, . . . , a„i) be an n x m matrix such that ^ for all i, and 
let b — {bi, . . . , &„i)^ e R™ be an TO X 1 constant vector. Define a closed 
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convex polyhedron by 

X = {x e M" I A^x < b}. 

Throughout the paper, K is assumed to be nonempty. Suppose that x G M" 
is a random vector distributed as an rt-dimensional standard normal distri- 
bution Nn{0,ln)- The primary motivation and hence one of the purposes 
of this study is to evaluate the n-dimensional normal probability 

P{K) = Pn{K) := Pt{x e K). (1) 

When m = n and A is non-singular, if is a cone referred to as the 
simple cone formed as the intersection of n half spaces located in the general 
positions (Barvinok, 2002^). In this paper, we do not require that the apex 
of cone lies at the origin. In the case where m < n and rank(A) — m, K is 
decomposed as K = Ki © ]R"~™ with Ki an m-dimensional simple cone, 
where "®" means the orthogonal direct sum. In this case the probability ([T]) 
is obtained as P„(A') = Pm{Ki). For such a simple cone K, Miwa, Hayter 
and Kuriki (2003^ proposed a recursive integration algorithm to evaluate 
the probability P{K) by generalizing the idea of Abrahamson (1964)."^ This 
algorithm is practically useful for dimensions n up to 20, and is available 
in the R library mvtnorm.'^ 

In this paper, using this integration algorithm as a building block, we 
demonstrate that the probability P{K) for any polyhedron K can be com- 
puted by means of the abstract tube ideas proposed by Naiman and Wynn 
(IQQ?).!^ This method provides a sophisticated version of the inclusion- 
exclusion identity, and divides the complimentary set K'^ into signed cones. 
By calculating the normal probability for each cone and summing them up, 
we obtain P(K) = 1 — P{K'^). However, as we shall see later, the resulting 
cones are not necessarily simple, and we need to introduce a perturbation 
into the linear inequality system. We focus on a particular type of pertur- 
bation, and refine the theory of abstract tubes for a perturbed inequality 
system in this setting. Some statements are simply taken from Naiman and 
Wynn (1997),'^ but some of them are novel. For example, we show that 
under our perturbation the perturbed system always defines an abstract 
tube instead of a "weak abstract tube" . 

Moreover, we show that the abstract tube for the perturbed system can 
be constructed by means of a standard linear programming technique. We 
propose that such an algorithm should be used, and we make some re- 
marks that should be incorporated into its implementation. The prototype 
R program is available from the authors. 
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The construction of the paper is as follows. In Section [21 the theory 
of the abstract tube for a perturbed system is investigated. In Section |31 
the algorithm and its implementation for constructing the abstract tube is 
discussed. In Section|31 we exhibit numerical examples. For Tukey's studen- 
tized range statistics, the abstract tubes are constructed in various settings, 
and their distribution functions are computed. 

2. Abstract tubes for a perturbed inequality system 

In this section we investigate the abstract tube associated with a closed 
convex polyhedron. Because of the reason stated later, we focus on the 
case where the inequality system defining the polyhedron is perturbed. We 
propose a lexicographic perturbation in an outer direction (see ([4])) and 
refine the theoretical results of Naiman and Wynn (1997^ in this setting. 
Define m half spaces 

H^^ {x eW \ ajx <b,}, i = l,...,m. 

The complement and the boundary of Hi are denoted by 

i7,f = {x e R" I ajx > b,} and dH, = {x G R" | ajx = b,}, 

respectively. Since P{K) = P{nYl^H,) = 1 - P{U'^^^Hf), we can evaluate 
^'(U™ i-fff ) instead of P(n™ li?,)- 

The family of the complements of the half spaces is denoted by 

H^{H^\t^l,...,m}. 

Let 

F,::^dH,r\K, i = l,...,m. 

If i^i 7^ 0, i^i is a face of K. Define a family of subsets of the indices 
{!,..., to} by 

^ = { J C {1, . . . , m} I J ^ 0, n.ejF, ^ 0}. 

Then, J- forms an (abstract) simplicial complex. {Hi^jFi \ J G J-} is the 
family of all faces (vertex, edge, facet) of i^T. Note that Hazj^Fi = Cii^j^Fi 
can happen even though Ji ^ J2. 
Let 

fl, iixe S, 

is[x) = < 

10, otherwise, 

be the indicator function of the set S. The following is Theorem 2 of Naiman 
and Wynn (1997) .1^1 



4 



Proposition 2.1. The pair ['H,T) is an abstract tube in the sense that 

= Y.^-l)\-'\-^ln^^,m{x) for all x G R", (2) 

where \ J\ is the cardinality of the set J . 

Assuming that x € R" is a standard Gaussian vector distributed as 
-^n(0,/n), and taking the expectation of ©, we obtain the formula 

l~P{K)=Y,{-lf'\-'P{n,e.,H'r). 

In the right side, if C]i^jHf forms an n-dimensional simple cone or the direct 
sum of a simple cone and a linear subspace, we can use the algorithm of 
Miwa et al. (2003|2!for evaluating the probability P{r\i,= jHf). However, this 
is not always the case, and for it to be the case the condition of "general 
position" is required. Before defining this notion, we explain two examples 
of abstract tubes that do not meet this condition. 

Example 2.1. Consider a polyhedral cone called a "pyramid" K C R'^ 
consisting of {xi, X2, x^) satisfying 

-Xi - X2 + X3 < 1 [1] 

—xi + X2 + x:i < 1 [2] 

+Xl + X2 + .X3 < 1 [3] 

+xi -X2+X3 <l [4] 

The boundaries of these half spaces are not in the general position in the 
sense that 4 hyperplanes share a point (0, 0, 1) in R'^. The simplicial complex 
is shown to be 

T = {1, 2, 3, 4, 12, 13, 14, 23, 24, 34, 123, 124, 134, 234, 1234}. (3) 

In the expression above "134" means {1,3,4} for example, and this con- 
vention is adopted throughout the paper. Applying this to the Proposition 
12.11 we can decompose 1 — P{K) into | — 15 terms. However, this decom- 
position is not appropriate for our purpose because the term 1234 stands 
for a non-simple cone (pyramid), and hence this decomposition does not 
simplify our problem. 
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Example 2.2. Consider the following system containing a redundant in- 
equality: 

+Xi — X2+ 0X3 < 
— Xi — X2 + 0X3 < 
Oxi — X2 + 0X3 < 

The third inequality is redundant. This 
dimensional, and 3 boundaries (lines) meet at 
cial complex is shown to be 

{1, 2, 3, 12, 13, 23, 123}. 

This decomposition is also unsuitable because the 3-dimensional term 123 
appears although this example is essentially 2-dimensional. 

To avoid unfavorable events such as have appeared in these examples we 
need the assumption that {dHi} is in the general position defined below. 

Definition 2.1. The set of m hyperplanes in R", for example, {dHi C 
M" I i = 1, . . . ,m}, is said to be in the general position when there does 
not exist J C {1, . . . , m} such that Cii^jdHi contains a max{n + 1 — | J|, 0} 
dimensional affine subspace. 

Remark 2.1. This definition is equivalent to that in Stanley (2007).'^ The 
definition by Naiman and Wynn (1997)^ is weaker than ours. In their defi- 
nition, when there does not exist J C {1, . . . ,m}, \ J\ = n + 1, Di^jdHi ^ 0, 
the set of hyperplanes is said to be in the general position. Example 12.21 is 
in the general position in their definition, whereas it is not in the general 
position in our definition. 

When a given {dHi} is not in the general position, we can apply an 
infinitesimal perturbation to rearrange the system into the general position. 
In this paper we restrict our attention to the lexicographic perturbation in 
an outer direction proposed below. 

For (bi) e M'" and e > 0, define 

b{e) = = (5, +e') = {bi+e\..., b„, + e'")^ 

(e* is e to the power i), and 

K{e) = {xeW \ A^x < b{s)}. (4) 

Define Hi{e), Fi{e) and Ti^e) similarly. Note that K = K{0), H, = H,{0), 
F, = F,(0) and J' = J"(0). 



[1] 

[2] 
[3] 

example is essentially 2- 
the origin in . The simpli- 
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The following two lemmas, Lemmas 12.11 and 12.21 are refinements of 
Lemma 1 of Naiman and Wynn (1997)® under this lexicographic perturba- 
tion and with the stronger definition of the general position. 

Lemma 2.1. For all sufficiently small £ > 0, 

(i) the family of hyperplanes {dHi{e) \ i = I, . . . , m} is in the general 
position, and 

(ii) F{e) does not depend on e. Write F{Q+) = F{e). J^(0+) is a 
simplicial complex. 

Proof, (i) Fix J C {1, . . . ,m}. Consider first the case \ J\ < n. Let Aj = 

{ai)i(zj be n X |J|, and let bj{e) — {bi{e))i^j be |J| x 1. Suppose that 
n,:gj9i/i(e) contains an n+1 — | J| dimensional affine subspace {x = Cd+xo \ 
d g M"+-'^^l''l}, where C = (ci, . . . , c„_|_i_|j|) is an n x {n + 1 — \ J\) matrix 
with linearly independent column vectors. Then Aj{Cd + xq) — bj{e) = 
for all c?, and hence AjC = and AjX^ — bj{e) = hold. Since rank(C) = 
n + 1 — I J|, we see that rank(^j) < | J| — 1. Hence Aj can be written as 
Aj — EF, where i? is n x r and F is r x | J| with r = rank(^j) < | J| — 1. 
Substituting this, we have 

= Ajxo - bj(e) = F^(S^xo) - bj{e) = {F'^ ,-bj{e)) 

Therefore, 

fj{e) := detj (_5^^^)t) (^^, -bj{e))^ 

= det(FF^) X - F'^ {FF'^Y^ F] bj{e) = 0. (5) 

Since rank(i^) = r, fj{e) is a polynomial in e of at least degree 2r(> 2). 
Then 

Sj = sup{(5 > I /,/(e) > 0, Ve e (0, S)} > 0, 

since the number of zeros of fj is finite. For all e G {0,Sj), r)i£jdHi{e) 
contains no (rt + 1 — | J|)-dimensional affine subspace. 

For the case | J| > n + 1, assume that xq e r\i^jdHi{e) exists. Since 
r = rank(v4j) < n < | J| — 1, we can follow the proof for \ J\ < n above. 
We define the polynomial fj{e) in and conclude that for all e € (0, 5j), 
Hi^jdHiie) = 0. 

Let (5* = minj Sj > 0. For all e G (0, (5*), {dHi[e) | z = 1, . . . , m} is in 
the general position. 
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(ii) The proof for (i) implies that for a given J C {1, . . . , m} the feasi- 
bihty or infeasibihty is unchanged for ah e e (0, (5,). □ 

Corollary 2.1. For each J e 

(i) {ai I i € J} is linearly independent, and 
(ii) the number oj elements \ J\ is at most rank(^). 

Proof, (i) Note first that Di^jdHi contains a point (i.e., 0-dimensional 
affine subspace), and because {dHt | i e J} is in the general position 
(Lemma |2.H (i)), it should be that max{n + 1 — 0} > 1 or | J| < n. 

Suppose that {ai \ i G J}, J G J-{0+), is linearly dependent. Let L = 
span{ai | i G J} and N = {x \ aj x = 0, i e J}. Then | J| > dim(L) + 1 = 
n — dim(A^) + 1 and hence dim(iV) > n + 1 — | J|. On the other hand, 
Di^jdHi ^ N + xq, where xq G Di^jdHi. This contradicts the fact that 
{dHi I z e J} is in the general position. 

(ii) The number of linearly independent vectors {ai \ i G J} must be at 
most dimspanjai} = rank(A). □ 

Lemma 2.2. Let J- and J-{0+) be simplicial complexes defined in Propo- 
sition [KT\ and Lemma \2.1[ respectively. Then J- I) J^(0+). Consequently, 
\T\ > |J-(0+)|. 

Proof. Suppose that J G J-'{0+). Wc prove that J G J-. From Lemma 
12.11 for all sufficiently small e > 0, it holds that Ej{e) n K{e) ^ 0, where 
Ej{e) = C^^(.J^H,{e) = {x | Ajx = b{e)} with Aj = (a^Ojgj (n x \ J\). 
Let Ej — Ej{Q). From Corollary 12. 1[ the column vectors of Aj are linearly 
independent, and hence Ej ^ 0. 

Assume that Ej DK — (d. Since both Ej and K are closed polyhedra in 
M", there is a hyperplane separating Ej from K, and the distances of the 
hyperplane from Ej and from K are greater than S > 0. By choosing e > 
such that 

■ Tr~ir < °' 

1=1, ...,m lloill 

this hyperplane also separates Ej{e) from K{e). Therefore, Ej{e)r\ K{£) = 
0, which contradicts the assumption that J G J-{0+). Hence, Ej D K ^ (d, 
from which J G F follows. □ 



The following theorem gives another version of the abstract tube as- 
sociated with a closed convex polyhedron K. This is an improvement of 
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Corollary 1 of Naiman and Wynn (1997)P who proved ^ for almost all x 
(i.e., a weak abstract tube). 

Theorem 2.1. The pair {'H,J-{0+)) is an abstract tube in the sense that 

lu™^^c(x)= (-l)''''"'ln.,.iff(a^) forallxeR'^. (6) 

JeJP(o+) 

Proof. As in Corollary 1 of Naiman and Wynn (1997) P it can be easily 
proved that ^ holds for x ^ U"^idHi. We show that © holds for every 
X G lJT=i9Hi. 

If a; e if n U"]^idHi, i.e., x is on the boundary of K, then x ^ for 
all i, and ^ holds as = 0. 

Suppose that x e K'^ DUl'LidHi. Recall that for sufficiently small e > 0, 

JeJ='(o+) 

holds. Since K is closed and e is sufficiently small, the left side of ([7]) 
becomes lu^^_H-,(e)c (x) = lu^^^Hf{x) = 1. 

When a; G a^a; — bi < bi{e), and hence lffj:(2;) = l//.(£)c(x) = 

0. When x ^ dHi, l//c(x) = 1//. (5)0(2;), because the distance between x 
and dHi is positive. Therefore, in the right side of d?]), l(-]-^jHi{e)''ix) = 

Summarizing the above, even when x G K'^ n U™ j^SiJ^, the expression 
([7]) with Hi{eY replaced by Hf holds, which is nothing but □ 

Remark 2.2. The abstract tube {H,F{Q+)) substantially improves the 
abstract tube {H,J-) because 

(i) the number of elements of T{0+) is not greater than that of T 
fLemma l2.2p . and 

(ii) DitzjHf, J G J-{0+), is a simple cone in M", or is the direct sum 
of a simple cone and a linear subspace fCorollarv l2.1l) . This is not 
always the case for J^. 

We present examples of abstract tubes for perturbed polyhedra below. 

Example 2.3. Applying a perturbation to the inequality system in Exam- 
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pie [211 

-Xi - X2 + X3 < 1 + [1] 

-xi+ X2 + X3 < 1 + [2] 
+xi + a;2 + ^3 < 1 + [3] 
+xi -X2 + xz<l + e^ [4] 

where e is an infinitesimal positive real number. The corresponding simpli- 
cial complex is shown to be 

J"(0+) = {1, 2, 3, 4, 12, 14, 23, 24, 34, 124, 234}, 

which is a proper subset of in ([3]). The number of terms |J^(0+)| = 11 is 
less than |J'(0+)| = 15. The maximum number of elements of J e J^(0+) 
is 3 (= n). The difference between F and J^(0+) is 

j-\ j-(0+) = {13, 123, 134, 1234}. 

The members of are characterized as the sets containing 13, that 

is, the non-existent edge in K{e). 

Example 2.4. Applying a perturbation to the inequality system in Exam- 
ple [H 

+xi - ^2 -f 0x3 < -h [1] 
-xi - X2 -)- 0x3 < -h [2] 
Oxi - X2 + 0x3 < + [3] 

The corresponding simplicial complex is shown to be 

^■(0+) = {1, 2, 3, 13, 23}. 

The largest |J|, J e J^(O-I-), is 2, as proved in (n) of Corollary [2Tl Note 
that a different perturbation obtained by altering the order of inequalities 

+xi - X2 + 0x3 < + [1] 
-xi - X2 -f 0x3 < -I- [2] 
Oxi - X2 0X3 < £^ [3] 

yields a different simplicial complex 

= {1, 2, 12}. 

This is an example where |J^(0+)| depends on the perturbation. 
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Taking the expectation of (jG]) with respect to the normally distributed 
random vector a:, we obtain 

i-piK)= i-if^-'pin.ejm)- 

As mentioned in Remark 12.21 ^iejHf appearing in the right side is of the 
class of cones whose probability can be calculated by the method of Miwa 
et al. (2003) P 



3. Construction of the abstract tube 

In order to construct the abstract tube (H, J^(0+)), we need to determine 
whether the system 

aj X ~ bi{e), i £ J, 
ajx < bi{e), i ^ J, 

has a solution for a given subset J C {!,..., m}, where e is an infinitesimal 
positive real number. For that purpose, linear programming (LP) techniques 
for degenerate systems are useful. 

First rewrite the problem in a linear programming format as follows: 

Problem 3.1. For a given subset J C {1, . . . , m} determine whether there 
exists a feasible solution x,y > (x,y £ K") and > (i = 1, . . . , m), 
'iJj ^ (j G J) such that 

Ui + aj {x — y) — bi{e) = 0, i ~ 1, . . . ,m, 
Vj - aj{x - y) +bj{e) ^0, j e J. 

To solve this problem we first define a tableau (matrix) corresponding 
to the above linear system. Let Aj — {aj)j^j (n x | J|), bj{e) — {bj{e))j^j 
(I J| X 1), and define an M x iV matrix 

[tij]MxN = 

where M = to + | J| + 1, iV = 2n + 1 and 1^ = (1, . . . , l)ixn- An algorithm 
for checking the feasibility of J is given as follows (Theorem 4-4-3 of Iri, 
197^). 



A' -A' -b{e) 
-Ajj Ajj b.j{e) 



Algorithm 3.1. 
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1. Choose an index i {—: ig) such that i < M — 1 and tiN > 0. // no 
such i exists, then the system is feasible. 

2. Among j 's such that j < N — 1 and ti^j < (if no such j ex- 
ists, then the system is not feasible) , choose an index j {—: jo) that 
maximizes tMj /\tio,j\- 

3. For i e {1, . . . , Af} \ {^o} and j G {1, . . . , iV} \ {jo}, let 

tij — tijgtif^j /tigjg, 

tija '■— tija/ti^jf^ , 

t'ioj '■~ ^tigj /tif^jg, 

tiojo ^/tiojo ■ 

4- Go to Step 1. 

In Step 1, tiN is a polynomial in e. ti^ > means that k{Q < k < n) exists 
such that tiN = Cfcs'' + Cfc+ie'^^^ + • • • + c„e" and Ck > 0. 

Using Algorithm 13. 1[ for every nonempty subset J C {!,..., m} we can 
confirm whether J is feasible (i.e., J G J-"(0+)). The following points should 
be incorporated into the implementation of the algorithm. 

Remark 3.1. Thanks to Corollary 12.11 the range of ,/ for checking the 
feasibility can be restricted to 

Jr,m = {J\JC{l,...,m},0<\J\< r}, (8) 

where r = rank(yl). Moreover, if Ji C J2 and Ji ^ J^(0+), then J2 ^ J"(0+). 
This reduces the range for checking the feasibility. 

Remark 3.2. In Algorithm 13 . 1 1 the polynomial ti^ = co + cif ^ + • • • + c„£" 
can be represented QbS cin 71 -|- 1 vector (cq i ci , . . . , Cn ). In this representation 
tiN > means (cq, ci, . . . , c„) > (0, . . . , 0) in the lexicographic order. This 
technique is known as the lexicographic method in linear programming.!^ 

Remark 3.3. In Algorithm l3. li the judgment that a value z is positive (z > 
0) should be implemented as z > eps, where eps is a small positive number 
comparable to the machine cpsilon. For example, the machine epsilon in R 
(32-bit) is 2^^^ = 2 x 10^^^, and in our prototype R program eps is set to 
10^14^ 
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Remark 3.4. In coding a computer program it is convenient to treat the 
candidate set Jr,m in ([51) as a totally ordered set. Note that 



We introduce a total order as follows: For J\ — {di, . . . (d\ < ■ ■ ■ < 

di^) and J2 = {ei,...,ez2} (ei < ••• < e;J in Jr.m, define Ji ^ J2 if 
h < h, or ^1 = Z2 and (di, . . . , d;J < (ei, . . . , e;^) (lexicographically). Then 
J — {di, di} {di <■■■< di) is the L-ih smallest element in Jr,m, where 



m — k 



L = i + u-i,™i+^ 

j — l k—dj-i-\-l 



i-j 



Here we let do = 0, and X^feLdJ-i+i = if dj-i + 1 > dj — 1. 

Conversely, for a given L such that 1 < L < \Jr.m\, the corresponding 
L-th smallest element J = {di, . . . , d/} (di < • • • < d;) in J7i-,m is obtained 
as follows: 

/ := max|/ > 1 | jJi-i,™] < i|, 
and for j — 1, . . . , / , iteratively, 

dj := maxjdj > dj_i + 1 | ^ f ) < ^ ~ |J;-i,m||- 



4. Tukey's studentized range 

The multiple comparisons procedure is a typical application where the mul- 
tidimensional normal probability in ([T]) is required. This is a statistical pro- 
cedure for combining several statistical tests. As an example, consider the 
comparison of k normal means. Suppose that, for z = 1, . . . , fc, is inde- 
pendently distributed according to a normal distribution with mean fii and 
variance ct|, and consider testing the equality of the means /i^, i = 1, . . . , fc. 
For X = {Xi, . . . ,Xk), Tukey's studentized range statistic is defined as 

T{X)^ max \^^1=£A_ (9) 

This is a combination of the test statistics \Xi — Xj\/ ^Jcrf + cr^ for testing 
/ij = iij. To determine critical values and power functions, we need the 
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probability of the event that the statistic (jQ]) takes a value less than or 
equal to a threshold c. 

The polyhedron K = {x E M.'' \ T{x) < c} is not proper in the sense 
that K contains a 1-dimcnsional linear subspace L ^ {x £ M.^ \ xi = ■ ■ ■ — 
Xk}- K is decomposed as K = Ki ® L, where Ki is a proper polyhedron. 
Fig. [T] is a picture of Ki when fc = 4 and the erf are all equal. When 
fc > 4 and the variances are equal (let af = 1), the facets of K are not in 
the general position because {Xi — X2)/V2 — c, {X2 — Xj^j^pl — c, and 
(X3 - X^)|^f2 = c imply {Xi - X^)/y/2 = c. One can see that 4 facets 
share a vertex in Fig. [T] 




Fig. 1. Studentized range (fc = 4) 



Table [T] shows the number of J^(0+) for fc = 2, . . . , 6 and constant vari- 
ances. Because Algorithm 13.11 depends on the order of the inequalities, the 
entries |J^(0-|-)| in Table [T] may vary when the order changes. However, nu- 
merical experiments suggest that the numbers |J^(0-f)| are independent of 
the order of inequalities. Moreover, even in cases where the variances of 
vary, |J^(0+)| seems to be unchanged. 



Table 1. Number of terms |J-"(0+)| in the ab- 
stract tubes 



fc 


2 


3 


4 


5 


6 


m = fc(fe — 1) 


2 


6 


12 


20 


30 


l-^(0+)| 


2 


12 


62 


320 


1682 



Note: m is the number of inequahties (facets). 
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In the statistical community there was a conjecture called the "Tukey- 
Kramer conjecture" which is stated as follows: Let Fk{ai, . . . ,af^;c) = 
Pr(r(X) < c), where X - iVfc(0, diag(crf )). Then for any k and c, 
Fk ((Ti , . . . , fTfc ; c) takes its minimum when the variances cr^ are all equal. 
This conjecture was affirmatively proved by Hayter (1984)!i2!Fig.[5]depicts 
the value F^lai, . . . , cts; c) with af = (10^)^*"^)/**, i = 1, . . . , 5, for the range 
s € [—5, 5]. The constant c is chosen such that i^5(l, . . . , 1; c) = 0.95. The 
minimum is certainly attained at s = 0, i.e., af = 1. 

Finally, note that although the polyhedron K depends on the threshold 
c, the shape of K is invariant for all c > 0. Once an abstract tube is obtained 
for some c > 0, this abstract tube is available for computing the distribution 
function Pt{T{X) < c) for any c in the null and non-null cases. 
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